Nonlocal transport mediated by nonlocal Hikami boxes: Condensation of evanescent 
quasiparticles injected into the superconducting gap 
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Evanescent quasiparticles entering a superconductor and propagating over distances larger than 
the coherence length give rise to intermediate states induced by double Andreev scattering on 
disorder. The resulting effective attractive interaction between evanescent quasiparticles is retarded 
at the extremely slow frequency of the applied bias voltage. The out-of-equilibrium mesoscopic 
superconductor with a fluctuating phase variable is compatible with a recent experiment [S. Russo 
et al., Phys. Rev. Lett. 95, 027002 (2005)]. Microscopic theory is discussed in the random phase 
approximation. 

PACS numbers: 74.45.+c,74.78.Na,74.78.Fk 
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I. INTRODUCTION 

The developments in nanotechnology allow to envi- 
sion the realization of a source of entangled pairs of 
electrons 1 ^, the electronic counterpart of a source of 
Einstein, Podolsky Rosen pairs of photons, with pos- 
sible applications to quantum information and to a 
test of non locality of quantum mechanics with elec- 
trons. Recent experiments on three terminal devices by 
Beckmann et al. 3 , Russo et al.— and more recently by 
Cadden-Zimansky and Chandrasekhar 5 realize an impor- 
tant step in this direction by probing crossed Andreev re- 
nectior^^ii .^^^^^^^^^^ ., in which Cooper 
pairs from a superconductor give rise to pairs of corre- 
lated electrons in two different ferromagnetic or normal 
electrodes. 

Andreev reflection^ 1 -, the mechanism by which an elec- 
tron from a normal electrode is reflected as a hole at a 
normal metal-superconductor (NS) interface while a pair 
is transmitted into the superconductor, takes place in a 
region of size £, the coherence length of the disordered 
superconductor. Crossed Andreev reflection (CAR) in 
a three terminal normal metal-insulator-supcrconductor- 
insulator-normal metal (N a ISIN&) trilayer corresponds 
to an Andreev process such that an electron from 
is transmitted as a hole in N a over a distance of order 
£, leaving a pair in the superconductor. Elastic cotun- 
neling (EC), the other competing channel, amounts to 
transporting an electron from N& to N a without chang- 
ing its spin. The two possibilities are shown schemat- 
ically on Figs. [1] and [H The nature of the domi- 
nant crossed transport channel can be controlled by 
the relative spin orientation of strongly polarized fer- 
romagnets in a ferromagnet-insulator-superconductor- 
insulator-ferromagnet (F a ISIFf,) structure^: EC (CAR) 
dominates with (anti)parallel spin orientations, as it can 
be seen from the spin of the electron or hole transmit- 
ted in electrode "a" (see Fig. [T]). On the other hand, 
EC dominates for normal metals with highly transparent 
interfaces^. Crossed transport dominated either by EC 
or by CAR is obtained in the three recent experiments 
mentioned above 3 ^ 4 ^. 
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FIG. 1: (Color online.) Schematic representation of the elec- 
trical circuit in a crossed transport experiment, and of the 
two possibilities for crossed transport in a N a ISINf, trilayer: 
transmission in the electron-electron channel (elastic cotun- 
neling, EC) and in the electron-hole channel, leaving a pair 
in the superconductor (crossed Andreev reflection, CAR), "e 
■f" corresponds to a spin-up electron and "h J." to a hole in 
the spin-down band. The CAR and EC transmission coef- 
ficient decay exponentially over the coherence length as the 
thickness of the superconductor is increased. 



More precisely, the crossed conduct anc o 22 ' 23 Ga,b{Vb) 
of a three terminal device measures the sensitivity of 
the current I a (V a , H) flowing through electrode "a" with 
voltage V a to the voltage V& applied on electrode "b" : 



Ga,b{Va,V b ) 



dIg(Va,V b ) 

dV b 



(1) 



(see the circuit on Fig. [5^). The voltage V a may be ar- 
bitrary but, according to available experiments 3 -^^, we 
consider that electrode "a" is grounded (V a = 0), as the 
superconductor. 

Theoretically, the crossed conductance defined by 
Eq. ((T|) vanishes to lowest order in the tunnel ampli- 
tudes for a NaISINf, trilayer because the electron and 
hole transmitted by EC and CAR have opposite charges, 
and because of a symmetry of crossed transport: CAR 
and EC have identical transmission coefficients^ ' 13 ! 14 (see 
Fig. [J) . This property does not hold for a single channel 



2 



INITIAL STATE FINAL STATE 




INITIAL STATE FINAL STATE 



FIG. 2: (Color online.) Schematic representation of the lowest 
order processes: crossed Andreev reflection (CAR, (a) and 
(b)) and elastic cotunneling (EC, (c) and (d)). Not only has 
the current through electrode "a" an opposite sign for CAR 
((a) and (b)) and EC ((c) and (d)) (as seen from the arrows 
on the figure), but the CAR current is exactly opposite to 
the EC current for the lowest order processes in the tunnel 
amplitudes (as seen from Ref.i). 



ballistic superconductor connected by single atom con- 
tacts to normal electrodes: in this case, the crossed con- 
ductance to lowest order in the tunnel amplitudes oscil- 
lates around zero with periodicity Xp (the Fermi wave- 
length) as the distance L between the contacts is varied. 
Oscillations in the Fermi phase factors average to zero 
in the realistic case of a multichannel contact 8 , or in the 
case of a diffusive superconductor 13 with non magnetic 
impurities. 

The theoretical prediction of a vanishing of the crossed 
conductance in a N a ISIN{, trilayer is however contra- 
dicted by the recent experiment by Russo et al£, that 
provides in addition evidence of a disappearance of the 
crossed signal in a magnetic field applied parallel to the 
layers. A characteristic energy scale friu c within the su- 
perconducting gap A is found experimentally 4 , at which 
the crossed signal changes sign from EC to CAR as the 
bias voltage energy increases above huj c , and eventually 
vanishes at higher energies. The experimental h\o c de- 
creases to zero as the superconductor thickness increases. 
One may conjecture that such experiments are character- 
istic of a geometry-dependent coupling of the supercon- 
ducting phase variable to crossed transport by evanescent 
waves. As we show, a natural possibility for such coupling 
(quantum interference effects in the superconductor) re- 
sults in an energy scale fku c smaller than the supercon- 
ducting gap for extended junctions. Another possibility 
(Coulomb interactions) was put forward by Levy Yeyati 
et al. in Ref. [Tfjl . in connection with the modes of the 
electromagnetic environment. 



II. QUANTUM INTERFERENCE EFFECT IN 
SUPERCONDUCTORS 

A. Transmission modes related to evanescent waves 

Weak localization enhances the return probability and 
the phase fluctuations of a normal metal^l. A supercon- 
ducting condensate is on the contrary delocalized and has 
a well defined phase. Superconductivity and weak local- 
ization are related to each other because the Cooper pairs 
of a superconductor can participate to the Cooperons of 
weak localization. Smith and Ambegaokar show indeed 
that the phase stiffness of a superconductor is reduced 
by weak localization 25 . Weak localization is however 
strongly modified in subgap transport through a super- 
conductor by the following property of evanescent wave 
functions. 

Electron tunneling through a normal metal involves 
the simultaneous forward propagation in time of an 
electron and the backward propagation in time of a 
hole, forming a diffuson. The terminology "transmission 
mode" is used here for describing the effect of disorder 
on subgap tunneling. For evanescent states, a hole prop- 
agating backward in time can be replaced in a transmis- 
sion mode by an electron propagating forward in time: 
the wave-function of a normal electron in a box of size £ 
being real valued, is equal to the time reversed particle 
wave-function in the absence of a magnetic field. Trans- 
mission modes of range £ in subgap transport through 
a superconductor can thus be made also either of a pair 
of electrons or of a pair of holes propagating forward in 
time2&. In a superconductor, a pair of electrons propa- 
gating forward in time can be obtained from a Cooper 
pair in the condensate. 

B. Electron- hole conversion in transmission modes 

As discussed in the preceding section, transmission 
modes in a disordered superconductor correspond to two 
quasiparticles scattering on the same sequence of impuri- 
ties. Transmission modes are described in a standard way 
on the basis of ladder diagrams (see for instance Smith 
and Ambegaokar— in the context of localization in a su- 
perconductor). As a direct consequence of Ref. HH, in- 
termediate states with electron-hole conversion from one 
impurity to the next such as on Fig. appear as the re- 
sult of disorder scattering, as opposed to usual Andreev 
reflection at a normal metal-superconductor interface 27 
being induced by spatial variations of the superconduct- 
ing pair potential. Intermediate states with electron-hole 
conversion in the disordered case are of course not due to 
anomalous electron-hole scattering on impurities, but to 
electron-hole conversion during propagation in between 
two impurities. 

For instance, transmission of an electron from elec- 
trode Nj, as an electron in electrode N a , and transmission 
of an electron in electrode Nf, as a hole in electrode N a 
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FIG. 3: (Color online.) Schematic representation of the ladder 
series capturing the influence of disorder on the processes on 
Fig. [T] Examples of Nambu labels ( "e" for electron and "h" 
for hole) are indicated on the figure, for elastic cotunneling 
with transmission of an electron ((a) and (c)) and for crossed 
Andreev reflection with transmission of a hole ((b) and (d)). 
"e j" corresponds to a spin-up electron and "h J." to a hole 
in the spin-down band. For clarity, the j and \ symbols of 
electrons and hole are not shown for the diffuson. 




FIG. 4: (Color online.) Schematic representation of trans- 
mission mode self-crossings in the presence of disorder. The 
circled area is the Hikami bo x 28 ' 29 , of dimension set by the 
elastic mean free path l e in a normal metal, and of dimension 
set by the coherence length £ for subgap transport. 



are shown on Figs. and b respectively. The extreme 
case of a sequence of electron- hole conversions due to dis- 
order, resulting in net transmission from electrode N& to 
electrode N a in the electron-electron channel is shown on 
Fig. [3J:. Another case of net transmission of an electron 
in electrode Nb as a hole in electrode N a by a sequence 
of normal scattering on disorder is shown on Fig. 
Electron-hole conversion in between two impurities on 
Fig. [3J; corresponds to a usual Andreev process, however 
"interrupted" at the time scale of the elastic scattering 
time, leading to intermediate virtual states with a charac- 
teristic energy much above the superconducting gap. We 
show below that this type of intermediate state charac- 
teristic energy scale for Andreev processes "internal" to 
the superconductor can be very much reduced by quan- 
tum interference effects and become observable. 



C. Self-crossings of transmission modes 

Self-crossings of transmission mode s 28 ' 29 are the build- 
ing blocks of weak localization in the normal cas o 24 i 30 
(see Fig. @J0- The distinguishing features of diffuson 
self-crossings in subgap transport by evanescent wave 
functions (see Figs. HJs and were already pointed 
in Ref. 0. The necessity of accounting for "advanced- 
advanced" or "retarded-retarded" transmission modes 
was already discussed by Altland and Zirnbauer— for an 
Andreev quantum dot. 

Weak localization- like loops such as on Figs.[4)D and [5)3 
contain, among all possibilities, intermediate states cor- 
responding to double Andreev processes in which an 
evanescent electron-like quasiparticle (an AR transmis- 



sion mode made of an advanced and a retarded Green's 
function for an electron propagating forward in time and 
a hole propagating backward in time) transforms by a 
first Andreev process at xi in an evanescent hole-like 
quasiparticle (another AR transmission mode for a hole 
propagating forward in time and an electron propagating 
backward in time) and the transmission of a pair (an AA 
transmission mode made of two advanced Green's func- 
tion for a pair of electrons propagating forward in time) , 
that recombine at X2 by another Andreev process after a 
propagation over |x2 — xi| ~ The intervention of the 
condensate in such loop processes is obvious from noting 
that the contribution of weak localization-like processes 
to the crossed conductance couples to a phase gradient 
in the superconductor—. These processes result in an at- 
tractive interaction of strength v per carrier injected in 
the gap, in units of the Fermi energy tp. This interaction 
due to quantum interference effects has all the features 
of a very slow phonon exchange: it is attractive, and re- 
tarded by a time IfmjeVb (because weak localization- like 
processes induces transitions only at energy hu = e14). 

As discussed previously, a transmission mode in the 
ladder approximation induces electron-hole conversion 
due to scattering on disorder, having the elastic scat- 
tering time as a characteristic time/inverse energy scale. 
Weak localization-like loops such as on Fig. QJ) corre- 
spond to the transport of pairs over a distance compara- 
ble to the coherence length £, and have thus the super- 
conducting gap A as a characteristic energy^. Similar 
double Andreev processes were discussed in Ref. [I?] in 
a NISIN three-terminal device, and, interestingly, were 
introduced previously by Jacobs and Kiimmel^ for the 
time evolution of wave-packets. As we show now, even 
smaller energy scales are obtained from many interact- 
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FIG. 5: (Color online.) Representation of (a) the electrical 
circuit for crossed transport through a N a SN(, trilayer, the 
electron e (top) and hole h (bottom) transmitted in electrode 
"b" by EC and CAR, and the "bare" transmission mode of 
wave-vector q in the ladder approximation. The junction has 
a dimension L y along the y axis perpendicular to the fig- 
ure. The aspect ratio is not to the scale of the experiment^ 
where L ~ 15 -f- 200 nm, and L x ,L y ~ 4, 8 /im. (b) is a sin- 
gle weak localization-like loop in a transmission mode self- 
crossing. The normal case is recovered for an AA transmission 
modes propagating locally over the elastic mean free path, (c) 
is a compact representation of (b) and of the associated tree. 
The four leaves correspond to the four transmission modes 
making a weak localization-like loop. AR (AA) transmission 
modes are represented by solid lines (with a cross), (d), (e) 
and (f) represent the random phase approximation (RPA) di- 
agrams with n = 2 loops and the associated trees. 

ing weak localization-like loops. The link between an 
effective dynamical model for a very slow attractive in- 
teraction in addition to the BCS pairing (Sec. and 
the resummation of many weak localization-like loops in 
the static case (Sec. lIV[) is shown schematically on Fig. El 



III. SLOW ATTRACTIVE INTERACTION IN A 
SUPERCONDUCTOR 

A. Energy scale within the gap 

Changing the phase by ±27r during a time inter- 
val 2tt/uj c requires the time 2tt/uj c to be smaller than 
2tt H/eVb, the retardation in the interactions that pre- 
vent the phase from fluctuating. The BCS case with 
a phase constant in time corresponds to retarded inter- 
actions with a Debye energy Hujd, orders of magnitude 
larger than the superconducting gap A. In the opposite 
limit eVb < Hu) c that we consider, the interaction v has to 



visit and to fix the phase of a number of state N(0)huj c 
in the normal metal to rotate the phase by 2tt within a 
delay 2tt/u> c , with N(0) — 2irLL x L y / X F ep the number 
of normal states per unit energy, and with tp the Fermi 
energy and Xp the Fermi wave-length (see L, L x and L y 
on Fig. [5]). By the uncertainty relation vepAt = h, the 
effective number of interactions in a time interval 2tt /uj c 
is Ni nt (uj c ) = 2ir/uj c At = 2-nvtp /huj c . The characteristic 
energy hw c is obtained by equating Ni nt (uj c ) = N(0)u> c : 

ef y LL x L y 

The decoupling between hu c and the states at energy 
larger than A is justified if fku c -C A, which we sup- 
pose in the following. The superconducting phase fluc- 
tuates by ±27r with a frequency oj c , but its average value 
over a time window larger than the retardation of in- 
teractions 2nh/eVb does not fluctuate. The Andreev 
current from the superconductor to electrode N a van- 
ishes for eVb < huj c because the superconductor chemical 
potential fluctuates between typical values ±Huj c . The 
crossed conductance is then dominated by EC, as it can 
be seen from second order perturbation theory in the 
tunnel amplitudes^ (see Fig. [TJ. 

B. Qualitative behavior of the crossed conductance 

We consider now the case eVb > hu> c . The proximity 
of the characteristic energy Hu> c induces pair correlations 
among evanescent quasiparticles injected in the super- 
conductor at bias voltage energies eVb > fiw c . CAR and 
EC processes are described schematically^: (b), (c) and 
(d) in the insert of Fig. [5] correspond to acting twice with 
the tunnel Hamiltonian, starting from the superconduc- 
tor in the BCS ground state for (b) and with a correlated 
pair in the initial state for (c) and (d). The crossed sig- 
nals due to the EC and CAR processes (c) and (d) cancel 
with each other because (c) involves an extra anticom- 
mutation of fermions compared to (d). The EC process 
(c) cancels with a CAR process (not shown on Fig. |5J) 
in which two electrons at energies — eVb and eV& in elec- 
trodes N a and N;, enter the superconductor as a pair. 
Fluctuations by one quasiparticle with spin-cr in the su- 
perconductor do not couple to the crossed signal because 
EC has an opposite sign for the two orientations of a. 
The remaining contribution to the crossed signal (not 
shown on Fig. |5J) is due to tunneling of the two electrons 
of a correlated pair from the superconductor to electrode 
N Q , leading to a negative differential crossed resistance, 
as for local Andreev reflection above the superconducting 
gap in a NIS junction 27 . From these arguments, we con- 
clude that a change of sign from a positive EC crossed 
resistance for eVb < to a negative Andreev reflec- 
tion crossed resistance for eVb > htu c occurs, and that 
the crossed signal disappears at higher energies, as in 
experiments^. 
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FIG. 6: (Color online.) Schematic representation the link between Sees. IIHI and IIVI A characteristic energy scale hui c is 
obtained from symmetry breaking between elastic cotunneling and crossed Andreev reflection (e ~ 1) in Eq. © and @. The 
energy scale hui c is identified to the characteristic frequency for phase fluctuations in the presence of a very slow attractive 
interaction, which, in turns, favors a symmetry breaking in favor of elastic cotunneling. 




(c) (d) 

FIG. 7: (Color online.) Schematic representation the lowest 
order processes for superconductor chemical potentials hui c 
larger than eVt ((a) and (b)), or smaller than eVb ((c) and 
(d)). The crossed Andreev reflection (CAR) processes (a) 
and (c) lead to opposite crossed currents for superconduc- 
tor chemical potentials ~ ±huj c while the elastic cotunneling 
(EC) processes (b) and (d) result in additive crossed currents. 
The crossed conductance with a superconductor chemical po- 
tential fluctuating between ~ ±huo c is thus dominated by EC. 



IV. RANDOM PHASE APPROXIMATION FOR 
WEAK LOCALIZATION-LIKE PROCESSES 

A. Method and algorithm 

1. Crossed conductance and resistance 

Now, we start from a static superconducting gap (see 
Fig. [5]) and determine from microscopic theory the ana- 
log of Eq. ||3J| for a dirty superconductor. We denote 
by T CA r(R, oj) and T BC (R, w) the CAR and EC dimen- 



sionless transmission coefficients of a superconductor at 
distance R and energy fko below the gap corresponding 
respectively to electron and hole transmission (see Figs.Q] 
and [2J described by transmission modes in the ladder ap- 
proximation such as on Fig. [3] 

We enforce a finite crossed current by breaking the 
symmetry Tbc(R)^) = Tcar(R-,w) according to 

T EC (R,u;) = T (R,w)(l + e) (3) 
Tcar(R, = T (R,w)(l-e). (4) 

The crossed conductance per conduction channel 

g a ,b(R, w) = -2 T 2 T c/l (R, u) (5) 

corresponds to processes of lowest order T 2 in an ex- 
pansion in the dimensionless interface transmission co- 
efficient < T < 1, and the charge transmission coef- 
ficient of the superconductor T c /j(R, to) = Tec^R, w) ~ 
Tc ar(R, accounts for the different carriers transmit- 
ted by EC and CAR in electrode "a". The resulting 
crossed resistance per channel is obtained as the inverse 
of the extra diagonal element of the inverse of the crossed 
conductance matrix: 



7? 


T^a,b 




Qa,a 


Ga,b 


-1 


_ Ttb,a 


T^b,b _ 




Qb,a 


Qb,b _ 





(6) 



with 

G ai , aj (Va,Vb) = -^(V a ,Vb) (7) 

generalizing Eq. (JTJ) , where the entries a* and dj in Eq. ([7]) 
correspond to the normal electrode labels "a" and "b". 
We deduce from Eq. ^ the off-diagonal matrix element 
of the crossed resistance matrix: 

K a ,b = —~ — r ^ a '\ — -p, — , (8) 

Va,aSJb,b — ya,byb,a 
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approximated as 

K a , b ~ -jfi^-, (9) 

y a , a yb,b 

because of the damping of crossed processes over the co- 
herence length £. The dependence on V a and Vi is not 
explicit in Eqs. © and We obtain from Eq. ([5]): 

Tl a , b (K, w) ~ ( T- 2 T C ,(R, w), (10) 

much larger than for our previous approach for higher 
order processes at the interface o 15 i 17 . 

2. Algorithm 

The excitations of the condensate, integrated out by 
the random phase approximation (RPA)21 for the trans- 
mission modes, lead to a new value of the crossed con- 
ductance corresponding to replacing e by e(R, w,e). We 
identify the characteristic energy hu c as the divergence in 
the energy dependence of e(R, cj, e), as for a gap edge sin- 
gularity. The symmetry breaking parameter e„(R, u>, e) 
with n localization loops is obtained by inverting 

f EC (n,e) = T (rc) + e n (e)T (11) 
TcAB,(n, e) = f (n)-e n (e)T , (12) 

with n the number of weak localization-like loops, and 
where the dependence on R and to is not explicit in 
Eqs. (HU) and (T2"]). We deduce the value of 

F n (q,w,e) = 1, (13) 

e 

where we changed variable from R to the wave-vector q 
by a Fourier transform. 

An algorithm generates the topologically inequivalent 
higher order RPA diagrams for the transmission mode 
with the tree structure shown on Figs. [5b, d, e, f. These 
diagrams maximizing the number of imbricated loops for 
a given number of branchings are the most relevant for 
describing intermediate states with a proliferation of An- 
dreev reflections internal to the superconductor for u> ~ 
lj c . In order to reduce the computation time, we first spe- 
cialize to a single channel by restricting to the transverse 
components (q' x , q' , q' z ) — (0,0,2Trn z /L) of the wave- 
vector q' = (q' x ,q y ,q' z ) = (2nn x /L x , 2im y /L y , 2irn z /L) 
in the weak localization- like loop, with n x , n y and n z 
three integers (see Fig. [S]). The scaling with the number 
of longitudinal components is restored afterwards. 

B. Results 

Now, we present the result of RPA resummation. The 
q-dependence of the transmission coefficients is, up to an 




J I I I I I L 



-3 -2 -1 1 2 3 4 5 
logicAo/l 

FIG. 8: (Color online.) Scaling plot of ^(A^e 2 ) as a function 
of A^e 2 for diagrams with n = 1 (•), n = 2 (a) and n = 3 (■) 

weak localization-like loops. The energy hu) and the dimen- 
sionless symmetry breaking parameter e vary over 3 orders 
of magnitude. The values of the elastic mean free path l e , 
such that 2 < < 14 (as compared to — 5 -f- 7 in the 
experiment^), vary from the dirty limit to the cross-over with 
the ballistic limit. For clarity, not all calculated points are 
shown on the figure, (a) shows one value 2-k/lu c of the super- 
conductor chemical potential, (b), (c) and (d) show EC and 
CAR processes for eVb > huj c with residual pair correlations 
among evanescent quasiparticles. 

overall rescaling, the same as for the bare transmission 
coefficient, because of the properties of the convolution 
of exponentials in the corresponding transmission coef- 
ficients in real space. We thus evaluate only the rescal- 
ing factor F„(0, u, e) that collapses on master curves (see 
Fig. [5]) when plotted as a function of the dimensionless 
parameter 

\^ = k F l(^j 2 (k F l e f\ (14) 

according to F„(0,w,e) = F n (X UJ e 2 ), with l e the elastic 
mean free path such that 2 < £/l e < 14, and with a — 
1.1±0.1. Summing the RPA diagrams leads to e, the new 
value of e modified by the response of the condensate: 

e>,e) = 6(l + F 0O (A LJ e 2 ))^ _\ 2 - (15) 



V. IDENTIFICATION OF THE ENERGY 
SCALES 

A. Characteristic energy of a disordered 
superconductor 

To take into account the large longitudinal dimensions 
of the junction, we note that convolutions of the expo- 
nential envelope of the transmission coefficients lead to 
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an enhancement of A^e 2 by a factor k F L x L y for extended 
interfaces. More precisely, combining diffusion modes for 
to a; and from x to L leads to convolutions of the form 




cxp (— x/£) exp (— (L — x)/£)dx = icxp (— L/£), 



( 16 ) 

where the one dimensional case with two evanescent dif- 
fusion modes is considered for simplicity. 

We deduce from Eq. (fT4| the characteristic energy huj c 
of a dirty superconductor: 

^£ = ±^ ■ 1 (17) 

e F |e| (Me) Q ^fk%LL x L y 

The parameter of interactions v in Eq. ([2]) is independent 
on A, as expected for a mechanism due to electron- hole 
conversion below the gap21; the lowest of the characteris- 
tic energy hw c and the voltage eVb, rather than ve F , is an 
indicator of the strength of weak localization-like quan- 
tum interference effect in subgap transport. The form of 
the transmission coefficients puts the constraint |e| < 1. 
The value e = 1 favoring EC with respect to CAR and 
minimizing fku c in Eq. (|17p for a static superconducting 
gap is compatible with the previous picture of a fluctu- 
ating chemical potential (see Fig. [6]), and with the large 
EC crossed signal measured experimentally at zero bias 
in Ref. H 

B. Comparison with experiments 

We use l e ~ 2nm from Ref. 0, e F ~ 5.3 eV, k F — 1 
A -1 for Nb. In experiments, hw c is limited by l v ~ 
0.1 -7- 0.2 fxm (as obtained from the inelastic electron- 
electron scattering time r e _ e ~ Ins 4 ), instead of the 
sample dimensions L x and L y in Eq. (|17|) . leading to 
fuv c ~ 21 10 fieV, 11-6 /xeV and 6 4- 3 fieV for L ~ 15, 
50, 200 nm respectively. Russo et alA find experimentally 
hoj c ~ 270, 50 /J.eV, for L ~ 15, 50nm, and fiu c below the 
resolution threshold for L ~ 200 nm. The theoretical 
approximation underestimates the experimental fku c in 
the regime L > £ (L = 15, 50 nm and £ = 10 -j- 15 nm 
in experiments^). The energy scale A is associated to 
objects of size < £ formed by the accommodated weak 



localization-like loops for L ~ £. The energy scale Sw c 
is thus enhanced by a cross-over to perturbative effects 
of weak localization- like loops, with huj c reaching ~ A 
for L ~ £, as in the absence of weak localization- like 
loop o 14 i 17 . The cross-over is not captured by Eq. (fl~7|) 
with L > £. 



VI. CONCLUSIONS 

In summary, intermediate states due to double An- 
dreev processes induced by a weak localization-like quan- 
tum interference effect in subgap transport are at the 
root of a very slow attractive interaction between evanes- 
cent quasiparticles. The interaction arises in an out-of- 
equilibrium situation if quasiparticles are forced to travel 
through the superconductor over distances exceeding the 
coherence length £, which does not apply to usual An- 
dreev reflection limited by £. The phase of the meso- 
scopic superconductor fluctuates with a characteristic 
time 2tt/lo c if the bias voltage energy eVJ, is smaller than 
Huj c , because of the large retardation h/eVb in double 
Andreev processes due to weak localization-like loops, 
as opposed to interactions retarded over the Debye fre- 
quency for a BCS superconductor. The model explains 
the following experimental facts: i) the existence of an 
energy scale Huj c within the superconducting gap, decay- 
ing to zero as the superconductor thickness L increases; 
ii) a change of sign from EC for eVb < hui c to CAR 
for eVfc > hu c , and a disappearance of the crossed sig- 
nal for Huj c larger than a few hio c ; and iii) the coupling 
of weak localization-like quantum interference effect to 
a magnetic field. A challenging issue is to account for 
the interplay between weak localization-like quantum in- 
terference effect discussed here and the coupling to the 
electromagnetic fielcU^. 
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